Faraday rotation in bilayer and trilayer graphene in the quantum Hall regime 
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Optical Hall conductivity, as directly related to Faraday rotation, is theoretically studied for bi- 
layer and trilayer graphene. In bilayer graphene, the trigonal warping of the band dispersion greatly 
affects the resonance structures in Faraday rotation not only in the low-energy region where small 
Dirac cones emerge, but also in the higher-energy parabolic bands as a sequence of satellite reso- 
nances. In ABA-stacked trilayer, the resonance spectrum is a superposition of effective monolayer 
and bilayer contributions with band gaps, while ABC trilayer exhibits a distinct spectrum pecu- 
liar to the cubic-dispersed bands with a strong trigonal warping, where the signals associated with 
low-energy Dirac cones should be directly observable owing to a large Lifshitz transition energy 
(- 10 meV). 



I. INTRODUCTION 

The observation of the anomalous quantization of 
(^xy = + ^) graphene quantum Hall system has 

established the existence of massless Dirac-quasiparticle 
in graphene^i^, and kicked off an increasing fascination 
with the massless Dirac physics of graphene. Optical 
properties have also been studied, where a point of inter- 
est is an unusual optical selection rule. Specifically, the 
optical Hall conductivity cFxyioo) is an interesting quan- 
tity to look at, since it is an ac-extension of the Hall 
conductivity which has a topological nature. 

The optical Hall conductivity has been theoreti- 
cally studied for graphene^—, graphit^ and topologi- 
cal insulator^. Experimentally, the optical Hall con- 
ductivity is measurable through Faraday rotation, since 
the Faraday rotation angle Qh is proportional to (Jxy{oo) 
in quantum Hall regime [where (no + ng) ^ \(yxx{^) i 
(Jxy{^)\/{ceQ)] as 
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where no{ns) is the refractive index of the air 
(substrate)-^, c the velocity of light, and Sq the dielectric 
constant of vacuum. Faraday rotation experimentally 
starts to be measured both for 2DEG^ and for graphene^^ 
in the quantum Hall regime. Due to the Dirac nature of 
quasiparticles in graphene, the optical Hall conductivity 
in graphene reflects the Dirac Landau level (LL) struc- 
ture {sn = sgIl{n)^/rlhuJc) and an unusual selection rule 
(|n| ^ |n| + 1), with Landau index n. 

Studies so far have focused on the optical Hall conduc- 
tivity in monolayer graphene quantum Hall system. In 
the physics of graphene, there are growing interests in bi- 
layer and trilayer graphene systems, since their electronic 
structures are distinct from that of monolayer graphene. 
For bilayer graphene the interlayer coupling between the 
two graphene sheets (71 in Fig(T]) changes the monolayer's 
Dirac cone into two parabolic bands touching at the K+ 
and K- points in the Brillouin zoneM'^ To be more pre- 



cise, if we take account of the second- neighbor interlayer 
hopping (73 in Fig(T]) a trigonal warping is induced on the 
band dispersion (Figdjb)). There, the parabolic bands 
touching at the Dirac points are reformed into four Dirac 
cones in the low-energy region. These are connected to 
the parabolic bands at a higher energy, at which a Lif- 
shitz transition in the topology of Fermi surface takes 
place. 

It is then an interesting question to ask how these will 
affect optical responses. While the optical longitudinal 
conductivity (Jxx{^) is discussed ini^J^, and experimen- 
tally utilized for determining the band structurei^ii^, we 
propose here that the optical Hall conductivity crxy{oo) 
should be a good probe for multi-layer graphene in the 
quantum Hall regime. Trilayer graphene, the next in the 
series of multilayer graphene, harbors another interest, 
since it comes with two different stacking orders, ABA 
and ABC. The low-energy dispersion of ABA trilayer 
consists of a monolayer- like Dirac cone and a bilayer-like 
parabolic bands, which are gapped due to the absence of 
inversion symmetry in the lattice structure ^i^Ti^^ On the 
other hand, ABC trilayer has a pair of cubic-dispersed 
bands with a larger trigonal warping effect than in bi- 
layer graphene ^ii^i^i22r— for which LL splitting due to 
Lifshitz transition is reported^. Thus the interest for 
trilayer is how the stacking dependence appears in the 
optical Hall conductivity. 

These have motivated us to study in this paper the 
optical (Hall as well as longitudinal) conductivities in 
bilayer and trilayer graphene systems. We shall show 
for bilayer graphene that the Lifshitz transition accom- 
panied by the trigonal warping greatly affects the res- 
onance structures in Faraday rotation not only on low- 
energy scale where Dirac cones emerges but also in the 
higher-energy range with parabolic bands as a sequence 
of satellite resonances, whose weight can become large 
in the vicinity of the Lifshitz transition, hence should be 
experimentally measurable. Thus, while it would be dif- 
ficult to directly observe the Lifshitz transition and low- 
energy Dirac cones due to the tiny energy scale of the 
trigonal warping (~ 1 meV) , the optical conductivities 
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FIG. 1: (a) Lattice structure of bilayer graphene with A/B 
sublattice carbons indicated in blue/red. The unit ceh con- 
sists of and Bi carbon atoms in the top layer and A2 and 
B2 carbon atoms in the bottom. 70: the nearest-neighbor 
(Ai ^ Bi,A2 ^ B2) hopping within each layer, 71: the 
vertical {Bi ^ A2) interlayer hopping, and 73: an oblique 
(Ai ^ B2) interlayer hopping. Two layers are separated by 
d = 0.334 nm. (b) Low-energy band dispersion of bilayer 
graphene with the trigonal warping effect with a Lifshitz tran- 
sition where the topology of Fermi surface changes from four 
Dirac cones to one, trigonally warped parabolic dispersion. 



should provide an experimentally accessible way to detect 
the trigonal warping as additional resonance structure. 

For trilayer graphene, on the other hand, we shall show 
that the optical conductivities are significantly affected 
by the difference in the stacking order. In ABA tri- 
layer, the resonance spectrum is a superposition of effec- 
tive monolayer and bilayer contributions with band gap, 
while ABC trilayer exhibits a distinct spectrum peculiar 
to the cubic-dispersed bands. In the latter, the trigonal 
warping effect is strong with a larger Lifshitz transition 
energy (~ 10 meV), so that we predict here the signals 
associated with low-energy Dirac cones should be directly 
observable. 



II. BILAYER GRAPHENE AND TRIGONAL 
WARPING 

A. Effective Hamiltonian 

AB-stacked bilayer graphene comprises two graphene 
sheets coupled by interlayer hoppings as shown in FigdJ 
where B sublattice in the top sheet is located just above 
A sublattice of the bottom sheet (Bernal stacking). A 
unit cell thus contains four carbon atoms, i.e., AI and 
BI in the top layer and A2 and B2 in the bottom. Hop- 
pings are the nearest-neighbor hopping within each layer 

(70) , the vertical interlayer hopping between BI and A2 

(71) . An oblique interlayer hopping between AI and B2 
(73) causes the trigonal warping of the band dispersion. 



70 = 3.2 eV, 71 = 0.39 eV, 73 = 0.32 eV. We shall adopt 
these values hereafter. There is another interlayer hop- 
ping parameters, 74 = 0.044 eV, connecting AI (BI) and 
A2 (B2), which introduces a small electron- hole asymme- 
try to the band structure^^i^ 

In a basis with components, V^ai, V^bi, '0A2, ^^^2, on the 
four sites, the Hamiltonian for the bilayer graphene is 
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where tt = ^TTa^ +i7ry, tt^ = ^^TTx — iTTy, tt = p + eA, with A 
being the vector potential arising from the applied mag- 
netic field, and the valley index = ±1 for K± points. 
Here v = \^ajo/{2h) is the band velocity for monolayer 
graphene, a ^ 0.246 nm the distance between the nearest 
A sites, and vs = V3ajs/{2h)^ V4 = ^/Za'^^/ {2%) veloci- 
ties related to 73, 74. 

The low-energy physics of bilayer graphene for a re- 
gion where the eigenenergy e is much smaller than the 
interlayer hopping 71, is captured by a 2 x 2 Hamilto- 
nian in a basis of A1/B2 carbon sites, with a perturba- 
tion with respect to 6/71, which is valid for energy scale 
\e\ <C 71 = 0.39eV, a^^^ 
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with an effective mass m = 7i/(2'u^). In the absence of 
magnetic fields, the first term on the right-hand side gives 
a pair of parabolic bands E = ±p^/(2m), while the sec- 
ond term coming from 73 causes the trigonal warping in 
the band dispersion, and the third term from V4 produces 
a weak electron-hole asymmetry by adding the band en- 
ergy 2vv4p'^/ji in both conduction and valence bands. 
In the low-energy region, two touching parabolae are re- 
formed into four Dirac cones as shown in FigHfb), where 
the Lifshitz transition (separation of the Fermi surface) 
occurs at 
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The Landau level spectrum in a uniform magnetic field 
B = rot A may be found with the relation (tt^tt^) = 
{V2h/£){a'^,a) for or (7r,7rt) = {V2h/£){a, a"^) for 
Here and a are raising and lowering operators, re- 
spectively, which operate on the Landau-level wave func- 
tion as a(j) n = V^^n-i,<^Vn = + l(/>n+i, while 
£ = y^h/ {eB) is the magnetic length. If we neglect the 

and V4 terms in EqnlSl the eigenenergies are given 



shujc \/n{n + 1), 



(5) 



The magnitudes of these are estimated in Refj 



33,34 



as 



where states are labeled by the Landau index n = 
— 1,0, 1,---, and the band index 5 = ± (defined for 
n > 1) labeling the conduction (s = +) and valence 
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FIG. 2: For bilayer graphene QHE system without the trigonal warping effect (vs = 0), we show (a) the optical longitudinal 
conductivity axx{eF,oj), and (b) optical Hall conductivity (cf, cj), grey-scale plotted against the Fermi energy er and the 
frequency uo. 



{s = —) bands. The cyclotron frequency ujc is the same 
as defined for a two-dimensional electron gas (2DEG), 
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Two zero-energy Landau levels (LLs) appear (n = 0, — 1) 
at each valley, while for large values of n the LLs tend 
to be proportional to n as in a 2DEG. The associated 
eigenstates for K-^ point are 



1 (n < 0) 
(n > 1) 



(7) 



where we defined 0^ = for n < 0. The wavefunction 
at K- is obtained by interchanging the first and second 
components in EqnH 

The ^3 term in Eqn. [3] causes a hybridization between 
Landau states ipn.s and ipn+?>m,s' •, while v/^ contributes 
to the energy shift of the Landau levels. The energies 
and eigenfunctions for high-energy Landau levels with 
|£^| > £^Lifshitz in > TiV2,/l {B > IT), are expressed 
in a lowest-order perturbation in vs and V4 as 
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In the low-energy region, \E\ < £^Lifshitz, on the other 
hand, the spectrum is reconstructed into monolayer- like 



Landau levels as ^at = sgn( A^) V SNhuJtrig for the cen- 
ter cone and sn = sgn{N)\^h(jOtrig for three off-center 
conesi^, where 
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The first excited level = 1 in this series appears in 
small enough magnetic fields such that hujtrig < T^Lifshitz, 
which amounts to B < 0.08T. 



B. Optical selection rules 

The optical longitudinal {(Jxx{^)) and Hall {axy{oj)) 
conductivities are evaluated from the Kubo formula as 
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where f{e) is the Fermi distribution, the energy of the 
eigenstate |a), j^^ = {a\ja\b) the matrix element of the 
current operator j = dH/dA^ and r] a small energy cutoff 
for a stability of the calculation which we set to 0.1 meV. 

In multilayer graphene the optical selection rule is of 
crucial interest, which just reflects the current matrices, 
Jx^jy^i appearing in the Kubo formula. Let us consider 
in the effective 2x2 Hamiltonian Eqnl3l for which the 
current matrices read 
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FIG. 3: (a) Low-energy band structure, and (b) Landau levels against magnetic field for bilayer graphene with the trigonal 
warping included, (c) The longitudinal axx{eF,oj) and (d) Hall axy{eF,oj) gre-scale plotted against the Fermi energy cf and 
the frequency u for a magnetic field B = IT (a vertical dashed line in (b)). (e) A diagram indicating allowed resonances in 
cTxy Vertical dotted lines indicate the Landau levels. Note that (n=-l) and (n=l, s=-l)" are different. 



The first term on the right-hand side of each equation 
is responsible for n ^ n ± 1 transitions. When we ex- 
press a = (n, s) and b = (n + 1, s') in the Kubo formula 
(Eqn ITT]) . the corresponding contribution is 

^.„,„+i_^.„+i,„ ^ _i^c„C„^, f-^{n + 1), (13) 

where we have dropped the band indeces s, s', since it is 
independent of their combination. 

In addition, the second term in the current operators 
in Eqn.[T2l and also the hybridization between n and n + 
3m in the eigenstates^^, give rise to additional transition 
series, 

n^n + l + 3m, (14) 
n^n + 2 + 3m, (15) 

with all the combinations of s and s' . A similar selection 
rule was previously found for graphitosis. In contrast, 
the third term including V4 in Eqn. [T2l contributes to the 



n ^ n ± 1 transistions as the first term, so that it does 
not cause any additional transitions but adds a small 
correction to Eqn. [131 

For the high-energy Landau levels with \E\ ^^Lifshitz, 
the contribution to (Jxy{oj) in the resonances of Eqns. 

(dH) and ([15]) becomes of the order of vl'^ and vl^^~^^\ 
respectively. In the lowest-order perturbation in vs^ the 
transitions are respectively allowed for n o n + 4 and 
n ^ n + 2, where the current matrices are respectively 
given by 

•n,n+4 -71+4,71 _ _ • f 
Jx Jy V 6 / ' 

.„,„+2.„+2,„^.(^|^^_ (16) 

Since the weights of these additional resonances have no 
dependence on the Landau index n while the n ++ n + 1 
resonance depends hnearly on n, the resonances arising 
from the trigonal warping should be relatively prominent 
in a low Fermi energy (ep) region. We also note that the 
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FIG. 4: Optical longitudinal [axx] dashed line) and Hall 
{axy'i solid line) conductivities plotted against the frequency 
u for a magnetic field B — IT and the Fermi energy ep — —20 
meV. 

factor jx^jy^ may have positive or negative sign depend- 
ing on the transition sequence. This property is specific 
to axy, since the corresponding factors in axx is \jx^\'^^ 
which naturally have no information about the sign of 
the resonances. 



C. Numerical results 

We start with calculating the optical Hall conductivity 
using Kubo formula Eqn(TT] for the bilayer system using 
the effective Hamiltonian (EqnISj) without the trigonal 
warping term. The result for the optical longitudinal and 
Hall conductivities in Figl2fa,b) shows that an intra-band 
transition (n, s) ^ (n + 1, 5) occurs around 

^^intra ~ ^^c- (17) 

This is natural, since LLs are almost equally spaced, un- 
like in the monolayer case where LL energy oc ^/n is not 
equally separated. On top of this, there are inter-band 
transitions across the band-touching point, with a se- 
lection rule (n, ±) o (n + 1,^). Thus the inter-band 
transition energy is 

^CJinter ^ 2|eF| (18) 

for large enough n. 

As stressed above, an important difference between axx 
and axy is that resonance factor j^j^^ is always positive 
for cFxx-, while it may have both signs for cFxy This puts 
a clear distinction for inter-band resonances of optical 
longitudinal and Hall conductivities. For cFxx^ resonances 
(n + 1, — ) (n, +) and (n, — ) ^ (n + 1, +) add up, 
and inter-band resonances occur over an entire region of 
Fermi energy between negative (n + 1, — ) and positive 
(n + 1, +) LL energies. By contrast, for (Jxy^ resonances, 
(n+1, — ) (n, +) and (n, — ) (^+1, +) have opposites 



signs in the resonance factor, so that they cancel with 
each other for a region of Fermi energy between (n, — ) 
and (n, +) LLs. So the inter-band resonance for cFxy is 
peaked at Cn+i,- < ep < Cn,- and en,+ < < en+i,+ 
with opposite signs. 

Now we examine the effect of the trigonal warping on 
the optical responses. We go back to the original 4x4 
Hamiltonian (Eqnl2]) including all the hopping terms in 
order to get quantitatively accurate results for all the 
energy regions. We have diagonalized the Hamiltonian 
with a basis set spanned by a finite number of Landau 
functions (/>n (0 < n < 100) to calculate the dynami- 
cal conductivities with the formula EqnlTTl In Figl3]we 
show a result for bilayer graphene in a magnetic field 
B = IT. A behavior of LLs with magnetic fields (Lan- 
dau fan diagram ) is indicated in FiglSl^b). At B = IT, 
the monolayer- like LLs in low-energy Dirac cones are not 
seen as the LL spacing already exceeds £^Lifshitz- 

Figures [3l^c,d) are (ixxi^F^^) and axyiep^^), respec- 
tively, plotted against the Fermi energy cf and the fre- 
quency uj. While the result becomes complicated, we can 
make an identification of the resonance structure in Fig. 
[3l^e), where we have extracted, from the Landau level 
spectrum (FiglSfb)) at IT, the expected positions of al- 
lowed resonances, as argued above. There we can actu- 
ally observe the satellite resonances of n ^ n + 2 and 
n ^ n + 4, which are absent in purely parabolic bilayer. 
It should be noted that, even though it is hard to di- 
rectly observe the Landau levels below the Lifshitz tran- 
sition due to the tiny £^Lifshitz, the satellite transitions 
in a region outside of £^Lifshitz should be observable as a 
manifestation of the trigonal warping effect. 

For axy{uj)^ the weight of the satellite peaks is con- 
stant as expected from the current matrices of Eqn. [161 
while that for the original peak (n o n + 1) depends 
on the Fermi energy through index n, as j^'^^^j^^^'^ — 
—in?i^/(2Pw?) :^ —leFji^m). So the relative weight of 
satellite to original peak is larger for smaller Landau in- 
dex n or smaller Fermi energy. 

To make the signs of the resonances clearer. Figure |4] 
plots the optical Hall conductivity Oxy against the fre- 
quency at a fixed Fermi energy ei? = 20 meV. We can 
see that, in addition to a large resonance n ^ n + 1 
around hujc = 3 meV, satellite resonances n ^ n + 2 
and n ^ n + 4 are seen. In the case of the longitudinal 
conductivity (Jxx the resonance appears as a symmetri- 
cal peak, while a resonance in (Jxy has an odd structure 
around a resonance frequency uJc Thus the sign in the 
resonance weight for axy determines the resonance shape, 
i.e., in accordance with the sign of Jx^jy^/i in Eqs. [T3land 
[T6l the resonance asymmetry of n ^ n + 2 is opposite to 
those of n ^ n + 1 and n n-^A. These features can be 
utilized for an experimental identification of resonances 
in the Faraday rotation, which is directly proportional 
to CFxyioo) with Eqnlll where the unit /h in Gxyioo) is 
related to 2a/ (no + n^) rad in Faraday rotation with the 
fine structure constant a 1/137. 

The resonance frequency for intra-band transition 
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FIG. 5: Lattice structure of trilayer graphene with (a 
stacking and (b) ABC stacking. 



ABA 



within the conduction band is larger than those within 
the valence band, which is a consequence of an electron- 
hole asymmetry due to 74 term. A deviation in the 
cyclotron mass for electron and hole bands prevents a 
complete cancellation between (n, — ) ^ (n + 1,+) and 
(n + 1,—) (^7+) transitions, which result in small 
interband transitions in a wide region of Fermi energy. 



III. TRILAYER GRAPHENE 

Now let us move on to optical responses of trilayer 
graphene QHE systems. Trilayer graphene occurs in two 
different stacking orders, i.e., ABA and ABC stacking 
orders as shown in FiglSl In both ABA and ABC trilay- 
er s, the spatial arrangement of the two successive layers 
is the similar to bilayer graphene, while the relation be- 
tween the first and third layers is different between two. 
The variation in the stacking order results in totally dif- 
ferent electronic structures, and here we show that this 
can be dramatically reflected in optical responses. 



A. ABA-stacked trilayer 

For ABA stacked trilayer, the effective Hamiltonian 
around Kj^/K- points is given by a 6 x 6 matrix (the 
dimension being 2 sublattices x 3 layers) a^ii^i^SiS^ 



Haba 



(19) 

where 70, 71 and 73 are the same hopping parameters 
as in bilayer. A' the on-site energy difference between 
the atoms with and without vertical bond 71, and 72(75) 
the next-nearest interlayer hoppings between Al and A3 
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{Bl and B3). We adopt the values for bulk graphite, 
72=-0.020eV, 75=0.038eV, A^=0.050eV^i^ 

With a unitary transformation^^—, this Hamiltonian 
is decomposed into two blocks as 



Ha 



BA 





Hb 



The first 2 by 2 block, 
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corresponds to a massive Dirac Hamiltonian with a shift 
in Fermi energy, while the other 4 by 4 block. 
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corresponds to a gapped bilayer Hamiltonian with an- 
other shift in energy. So the low-energy physics of ABA 
stacked trilayer graphene is effectively described as a 
superposition of gapped monolayer and gapped bilayer 
band contributions as seen in the low-energy band struc- 
ture in Figini^a). 

In FigEl^b), LLs are plotted with magnetic field B, 
which are labeled with M (B) for monolayer (bilayer) 
blocks. Landau index n, and band index ± (for n > 1). 
For moderate magnetic fields 5 ~ 1 T, the Landau level 
spacing oc \fB for monolayer is much larger than that for 
bilayer oc B. Since bands are gapped, zero-energy LLs of 
monolayer and bilayer appear at the bottom (top) of the 
conduction (valence) bands for Kj^ valley. 

We now look at the result for the optical longitudi- 
nal (Jxxi^F^^) and Hall conductivity (Jxy{^F-,^) plotted 
against the Fermi energy ep and frequency uj for ABA 
stacked trilayer graphene QHE system in Figsl6l^c,d). We 
can discern contributions from monolayer- like Dirac LLs 
and from bilayer LLs, both of which exhibit intra-band 
and inter-band transitions. Since Dirac cone is massive 
and MO LL is situated at the bottom of conduction band 
for Kj^ valley and the top of valence band for K-^ so 
that MO M1+ resonance occurs at a lower energy 

than Ml > MO for and vice versa for K- valley. 

A cancellation of resonances in cjxy^ due to opposite signs 

in current matrices, occurs between Ml > MO for Kj^ 

and MO M1+ for K- for a region of Fermi energy be- 
tween MO(i^+) and MO(i^_), while this is not the case 
with (J For bilayer contributions satellites appear due 
to the trigonal warping effect as (n, s) ^ (n + 1 + 3m, s') 
and (n, 5) ^ (n + 2 + 3m, s') as in optical responses for 
bilayer graphene (Fig|3j). 

So the message here is the ac optical responses in the 
ABA-stacked trilayer graphene accommodates a curious 
mixture of contributions from an effective massive mono- 
layer and from an effective gapped bilayer with the trig- 
onal warping effect. 
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FIG. 6: (a) Low-energy band structure, and (b) Landau levels against magnetic field for ABA-stacked trilayer graphene. (c) 
Longitudinal (Txx{eF^oj) and (d) Hall (Txy{eF^oj) plotted against the Fermi energy cf and the frequency uj for a magnetic field 
B = IT (a dashed line in (b)). (e) A diagram indicating allowed resonances in Gxy 



B. ABC-stacked trilayer 

If we turn to ABC stacked bilayer graphene, the ef- 
fective Hamiltonian around K point is a 6 x 6 matrix 
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We can derive a low-energy effective Hamiltonian as a 
2x2 matrix with basis for Al and B3, where we eliminate 
the states coupled by 71. As in the case of bilayer, a 
perturbation in s/^i gives the effective Hamiltonian for 
ABC trilayer graphene as^. 
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where tt^ = {tt^tt + 7r7r''")/2, and we neglected V4 term 
which gives the electron-hole asymmetry in a similar 
manner to bilayer graphene. In the absence of mag- 
netic fields, the first term gives a pair of cubic-dispersed 
bands touching at zero energy, while the second term in- 
volving 72 and vs causes a trigonal warping in the band 
dispersion. In a low-energy region, the cubic bands are 
reformed into three Dirac cones at off center momenta lo- 
cated in 120° symmetry around K± point. The Lifshitz 
transition occurs at 



72 

^'Lifshitz - Y 10 meV, 



(25) 



which is an order of magnitude greater than in bilayer's. 

If we first neglect 72 and vs and consider the cubic part 
alone in Eqn. [24l LLs are^^s^ 



Sn,s = s?iujABC\/n(n + l)(n + 2), (26) 
where n > —2 is the Landau index, s = ± (only for 
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n > 1) the band index, and 



ABC 



£ 



7i 



(27) 



A peculiar magnetic field dependence of cyclotron energy 
oc 5 2 for ABC trilayer graphene implies a smaller LL 
spacing compared to the monolayer's LL ex ^2 and bi- 
layer's LL oc B for weak magnetic fields. The associated 
wavefunction at point is 



(Pn-1 



(28) 



where Cn is defined in Eqn. [71 The wavefunction at K- 
is obtained by interchanging the first and second compo- 
nents in Eqn. [28l 

Similar to bilayer, the trigonal warping effect due to 72 
and vs hybridizes ?/^n,s with '0n+3m,s'- In the low-energy 
region |£^| < £^Lifshitz, the spectrum is reconstructed into 
the monolayer- like Landau levels from small Dirac cones 
as a series En = sgn{N) y/Nhuj trig^ where 



trig 



= 3 



72 


2/3 " 


271 





4^3 

3v 



271 
72 



1/3" 



1/2 




vV2heB 



(29) 



Each monolayer- like level is three-fold degenerate reflect- 
ing the number of small Dirac cones. For larger magnetic 
fields the degeneracy of the 7V-th level is lifted when en 
exceeds £^Lifshitz, and the three levels eventually connect 
to n = 3A^, 3A^ — 1 and 3N — 2 in the cubic band in 
Eqn. [26l The first excited level = 1 can appear when 
hujtrig < £^Lifshitz or 5 <^ IT, which is much larger than 
in bilayer graphene. 

Now we go back to the original 6x6 Hamiltonian in- 
cluding all the hopping terms (Eqnl23|) to discuss op- 
tical conductivities calculated with diagonalization and 
Kubo formula. In a low-energy region 10 meV), 
Dirac cones appear due to trigonal warping effects as seen 
in the band dispersion (Fig|71^a)). In the LL spectrum 
in Figd^b), a Lifshitz transition is clearly identified at 
^Lifshitz = lOmeV, which separates Dirac LLs and ABC 
cubic LLs. 

Figures [7Jc,d) show the optical longitudinal (Jxx{^f^^) 
and Hall conductivities (Jxyi^F^^) plotted against the 
Fermi energy cf and frequency uj for ABC stacked tri- 
layer graphene QHE system with B = IT. At this 
particular magnetic field, there are three Landau levels 
N = —1,0, 1 inside each Dirac cone, where the A" = ±1 
levels are close to the Lifshitz transision. The higher lev- 
els are outside of the Dirac cones, and can be regarded 
to belong to the the cubic dispersion. Accordingly we 
see cyclotron resonances between 1—^0 and 0^1+ 
in Ef < £^Lifshitz, while outside we see transitions be- 
tween cubic LLs (X 5 2 with much smaller level spacings 



than for Dirac LLs (x yB. In addition to resonances 
n ^ n + 1, the trigonal warping again gives rise to satel- 
lite transitions n ^ n + 4 and n ^ n + 2 with opposite 
resonance weights. There also emerge small resonances 
between Dirac n = LL and cubic LLs across the Lifshitz 
transition energy. Due to the ABC trilayer LLs arising 
from the cubic dispersion (Eqnl26|) the intra-band tran- 
sition energies show behaviors cx n^/^ with Landau index 
n, which is different from monolayer (oc n~^/^) or bilayer 
(constant), while the inter-band transition energies are 
qualitatively similar to ~ 26^. 

Thus an message in the case of ABC trilayer graphene 
is that the trigonal warping effect is an order of mag- 
nitude enhanced than in bilayer graphene, so it will be 
experimentally more feasible to access the effects induced 
by the trigonal warping. 



IV. SUMMARY 

We have studied the optical Hall conductivity for bi- 
layer and trilayer graphene. In bilayer graphene the trig- 
onal warping effect causes a Lifshitz transition with an 
emergence of four Dirac cones in a small energy scale 
1 meV. The trigonal warping is shown here to af- 
fect optical transitions as additional transitions arising 
from mixing of LLs whose Landau indices differ by 3 
and the opposite sign of resonance weights for transi- 
tions n ^ n + 2. For trilayer graphene, which comes 
with two stacking orders, we have shown that the opti- 
cal conductivity of ABA in low-energy region comprises 
contributionf from of monolayer- like massive Dirac and 
gapped bilayer with an electron-hole asymmetry, while 
ABC stacking has an order of magnitude larger trigonal 
warping effect, exhibiting resonances of Dirac LLs inside 
the Lifshitz transition (~ 10 meV), and a sequence of res- 
onances for the cubic bands with cyclotron energy oc B^ , 
This is predicted to be observable. 

There are several reports that the electron-electron in- 
teraction opens an band gap of a few meV at charge neu- 
trality point in bilayer and trilayer graphene^i^. If we 
consider the effect of the many-body induced gap in the 
effective models phenomenologically, this would cause a 
valley splitting of n = LL and energy shift of resonances 
associated with low- lying levels, while the transitions as- 
sociated to high LLs away from the Dirac points, includ- 
ing those of satellite transitions, should not be much af- 
fected. We leave the calculation including those effects 
for future works. 

There are various other problems. One is effects of 
disorder, which in this paper is treated as a phenomeno- 
logical broadening of the energies, assumed to be small. 
However, the localization physics is important for un- 
derstanding the quantum Hall physics, so that it should 
be interesting to treat disorder more accurately with 
methods such as diagonalization technique as done in 
Ref4^i^, and consider the localization effects on multi- 
layer graphene in the quantum Hall regime. It would be 



9 




Fermi energy (meV) Fermi energy (meV) Fermi energy (meV) 



FIG. 7: (a) Low-energy band structure, and (b) Landau levels against magnetic field for ABC-stacked trilayer graphene. (c) 
Longitudinal crxcc(eF, cj) and (d) Hall crxy{eF,(^) plotted against the Fermi energy cf and the frequency u for a magnetic field 
B = IT {a dashed line in (b)). (e) A diagram indicating allowed resonances in axy. 



another interesting problem to ask how various types of 
disorder, chiral-symmetry preserving and non-preserving, 
affect the optical responses in the multilayer graphene, 
since it is known that the presence or absence of chiral 
symmetry greatly affects the behavior of zero energy LLs 
for the bilayer graphene^. 
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